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Abstract 

It has been argued recently that the so-called nuclear modification 
factor (Raa) is an observable useful for identifying the quark-gluon 
plasma. We discuss the effect of Fermi motion in nuclei on Raa at 
CERN SPS and BNL RHIC energies. Contrary to the simple intuition, 
rather large effects are found for CERN SPS. The Fermi motion in 
nuclei contributes significantly to the Cronin effect. The effect found 
is qualitatively similar to the one observed experimentally at CERN 
energies and similar to the one obtained in the models of multiple 
scattering of initial partons. We predict different size of the effect for 
different types of hadrons, especially at low energies. 



1 Introduction 

It was realized only recently that the nuclear modification factor defined as 



R h AA (y,p t ;W) 



N, 



may be a very useful quantity which could shed more light on the dynamics 
of high-energy nuclear collisions PU El IH 03 E|- in the above definition 
y, pt and W are hadron rapidity, transverse momentum, and nominal NN 
center-of-mass energy, respectively. Nbi n is the number of binary nucleon- 
nucleon collisions usually calculated in the Glauber formalism. This brings 
in a small model dependence in the above definition. The index h indicates 
the measured hadron. 

Quite different values of nuclear modification factors have been observed 
at CERN SPS and at the Brookhaven RHIC 0111111111112!. It is speculated 
that the inequality 

R h AA (RHIC) < R h AA (SPS) (2) 



is due to the creation at RHIC of a dense nuclear medium, a new state of mat- 
ter, perhaps the long seeked quark-gluon plasma [U EE [H| . Various nuclear 
effects like shadowing, parton multiple scattering or final state quenching 
|3 El El modify the dependence of Raa on the transverse momentum and 
collision energy in quite different ways. In Ref.|2] the data from CERN SPS 
were interpreted within the picture of Glauber multiple scattering of initial 
partons (see (7j ) . It is not easy to make a priori predictions based on this ef- 
fect as one has to model parton-nucleon interaction. The proposed schematic 
models are flexible enough so that a rough description of the data is possible. 

In a partonic approach improved for including internal transverse mo- 
menta PUE], the average transverse momentum of partons in nuclear colli- 
sions is increased compared to the one in the elementary collisions 

(p t V + 6, ••.)• (3) 

The correction term used in the phenomenological analyses is typically a 
fraction of GeV, i.e. of the order of the nuclear Fermi momentum. 

At intermediate energies the inclusion of momentum distributions in nu- 
clei quantifying the nucleon motion is crucial for the production of kaons 
[TH| Uni HI] , eta mesons [IE] and antiprotons [THj in p+A collisions. It is well 
known that at subthreshold and near-to-threshold collisions the Fermi mo- 
tion leads to a huge enhancement of hadron production cross sections based 
on the first collision model. At the lower energies the inclusion of the off-shell 
effects incorporated through spectral functions (see for example [2lH I21[ l2~2*j) 
can be numerically important. 

At large energies one is tempted to completely ignore the effect of the 
Fermi motion of nucleons in the colliding nuclei. According to our knowl- 
edge, however, no quantitative estimates have been presented so far in the 
literature. 

In the present Letter we wish to show that the Fermi motion effect on 
R-AAipt) is not negligible even at high energies (W > 15 GeV) and must 
be included into consideration, e.g. in the context of the discussion on the 
creation of quark-gluon plasma. 

2 Folding model 

In order to demonstrate our main point, in the following we shall neglect all 
nuclear effects except the nucleon Fermi motion. We assume that the collision 
of nuclei A\(Z\, N\) and A 2 (Z 2 , N 2 ) can be considered as the collision of two 
ensambles of nucleons. In general, the production cross section of a hadron 
h depends on the type of the collision (a, (3 = pp, pn, np or nn). We assume 



2 



that at these high energies all off-shell effects can be neglected. In this 
somewhat simplified picture the cross section per elementary collision 1 is 

rf^AA-^h p iu 

a,/3=p,n ^ 

dy*dY t 8{y - y) 5\p t - p t ) W*) . 

In the above equation y*, Pf, W* denote rapidity, transverse momentum and 
energy in the nucleon-nucleon center-of-mass system, respectively. The quan- 
tities y and % (overall CM system) are obtained from pi, p 2 , y*, etc., via 
appropriate Lorentz transformations and rotations of the frame of reference. 
In Eq.Q pi is the nucleon momentum in the rest frame of the nucleus A\, 
while p 2 is the nucleon momentum in the rest frame of the nucleus A 2 . Eq.(J3J) 
is valid provided the standard normalization of momentum distributions is 
used 

fa/i(Pi) d A pi = 1 . (5) 



In the master formula (|%|). N a> p is the effective number of the collisions be- 
tween nucleons a and (3. Obviously N tot = ^ a „ N a> p. In the most naive 
approach N a , p /N tot = {N a / x Np/ 2 )/{A 1 A 2 ). 

In our intentionally simplified approach we consider the nucleus-nucleus 
collision as a collision of two bunches of nucleons. For simplicity we make 
no distinction between proton-proton, proton-neutron, neutron-proton and 
neutron- neutron collisions 2 . In addition, we shall make simplifications of 
kinematics relevant at high energies, when collective velocities of both nuclei 
are much larger than the Fermi velocity, i.e. va 1 /cm, v a 2 /cm 3> Vp. In this 
case Vi/ NN ,v 2 / NN || z \\ va 1 /cm-,v A2 / C m and y* « y. Then, the invariant 
distribution (per elementary collision) of hadron h in nuclear collisions is 
given as 

-{y, Pt] W) « / d 3 Pl d 3 p 2 f 1 (p 1 )f 2 ( P 2) -t-ptt^^ W) • (6) 



dyd 2 p t ' ' J dy*d 2 p* 

In the formula above, W* denotes the energy available in the elementary 
nucleon-nucleon collision. In order to calculate W* one needs to make the 
appropriate Lorentz transformations from four momenta in the rest frame 



1 This corresponds to dy ^ Pt in the definition of the nuclear modification factor (JJJ. 
2 This means that our present approach is valid only for the cases when: <J pp ^h z 
7 pn ^h — &np—>h = &nn— *h- This is roughly true for example for h= 7r°, ir + + 7r~, r\. 
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of nuclei [Ei/A^Pi/A^ to four-momenta in the overall center of mass of the 
colliding nuclei {Ei/cMiPi/CM)- Then the nucleon-nucleon cm. energy is 
calculated as 



W* = y {EyCM + E 2 /Cm) - (Pl/CM + P2/CM ) • (7) 

In the equation above the on-shell relation E^cm = \Jp1/cm + m N (^=1 ,2) 
is used. Due to plateau at midrapidities, at y « it seems reasonable to 
make an extra approximation y* = y. 

Although the master formula (JBJ) is very simple, its practical application 
requires knowledge of nucleon momentum distributions (fi and f'2) in nuclei 
as well as the rapidity, transverse momentum and energy dependence of the 
elementary production mechanism. It is needless to say that realistic calcu- 
lation of both nucleon momentum distribution in nuclei and differential cross 
section for elementary collisions is rather difficult. While the first requires 
nuclear many-body approach to the nucleus structure, the latter involves 
nonperturbative QCD effects which are not fully understood. The exact nu- 
clear calculation based on modern nucleon-nucleon interactions are available 
only for the deuteron, 3 if , 3 He, and quite recently - the four-nucleon system 
4 He [21]. Similarly, our understanding of the soft sector (small transverse 
momenta) of particle production within QCD is still not complete |23j . 



3 Estimate of the effect 

Before we go to the presentation of predictions based on Eq.©, we shall 
shortly describe the ingredients which enter the formula. In the Fermi gas 
model 

f(p) = C FG 6( PF -p), (8) 

where p F is Fermi momentum and 9 is the standard step function. The nor- 
malization constant Cfg can be obtained from the normalization condition 
(JSJ). In numerical calculations we take p F = 1.37 fm -1 . In the present paper 
we shall also consider a simple Gaussian parametrization of the momentum 
distribution of nucleons in 12 C nucleus from Ref.|19j. While the first seems 
more adequate for heavy nuclei, the latter may be representative for light 
nuclei. 

The elementary cross sections for the production of tt^, 7r°, in the 
proton-proton or proton-antiproton collisions are know experimentally in 
the broad range of energies. As mentioned above, in the present analysis 
we shall limit ourselves to the production of it only. For our purpose a 
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simple parametrization of the available world data would be useful. Such 
a parametrization exists, however, only for charged particles in the limited 
range of energies (2^]. In the parton picture, assuming isospin symmetry, the 
invariant cross section for the pp — > tt°X reaction can be expressed as 

do pp ^° , „ rx 1 / da pp ^ w+ , TTT , da pp ^ w ~ , „A 

■W ( "' ; = 2 (u^T ( "< ; + ( " ,; H/) J • (9) 

The effects which violate the above relation are expected to be sufficiently 
small. In the present calculation we use a parametrization of the ISR data 
from Ref-PS!- 

In general, there are three effects contained in the original convolution 
formula (j3J) which may result in Fermi motion modifying of R\a- The first 
effect is connected to the fact that p* t ^ p t . The second effect is related to the 
fact that y* ^ y. Because the general Lorentz transformation from the NN 
system to the overall nucleus-nucleus CM system includes also rotations of 
the frame axes, the formulae are rather cumbersome and not easy to reverse. 
Therefore these two effects cannot be included exactly in a simple way. The 
Monte Carlo method seems more suitable for this purpose. A Monte Carlo 
calculation, to be discussed in details elsewhere, shows that while the first 
effect is rather small, the second effect is almost negligible. 

The main (third) effect is connected to the fact that the true cm. en- 
ergy in the elementary nucleon-nucleon collision is not equal to the nominal 
nucleon-nucleon energy W* ^ W, even when averaged over Fermi distri- 
butions. As an example, in Fig[T] we show two distributions of the true 
nucleon-nucleon center of mass energy f(W*) 3 for two nominal energies W 
= 20 GeV (left panel) and W = 50 GeV (right panel). We observe a broad- 
ening of W* distribution with raising nominal energy, although the relative 
dispersion AW*/W stays roughly constant. The invariant elementary cross 
sections for particle production rise across the region over which f(W*) is 
spread (see FigEJ) leading to the value of (W*) c 4 larger than W. This effect 
is stronger for larger p* t and smaller energies W*. 

In FigOlwe show R\a as a function of the ir° transverse momentum for 
the 7T° rapidity y=0 and two different Fermi momentum distributions: 

(a) Fermi gas model distribution given by Eq.(|5|) (dashed) and 

(b) the Gaussian parametrization from Ref.^Hj (solid). 

The experimental data points were obtained by dividing the nuclear data 

3 J f{W*) dW* = 1. 

4 The index a indicates averaging of W* with the energy-dependent NN ~ * h cross 
section. 
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from Ref. [S] by the cross section for the elementary collision obtained with 
the help of Eq. © and the parametrization from [23] . 

The Fermi motion effect calculated from our master formula (jBJ) is strongly 
energy dependent, which is demonstrated in Fig0] Generally, the larger 
energy the smaller the effect. Surprisingly, the effect of the Fermi motion 
stays present even up to W = 50 GeV. 

In the present analysis we have estimated the effect of Fermi motion 
assuming that nucleons can be treated as free on-mass-shell particles. Under 
this assumption the light-cone fraction carried by a nucleon, £, defined so 
that 0< £ <A, is £ = (^/iti 2 n + p 2 + p 3 )/(iriA/ A). As a result the average 
< £ > is larger than 1 and the momentum sum rule is violated - nucleons 
carry momenta larger than the total nucleus momentum. Effectively the 
energy in NN collisions is overestimated. This is connected with the effect of 
nucleon binding in nuclei. 

The problem can be even better seen in the nucleus rest frame, where the 
total nucleon-gas energy is larger than the nucleus mass 

A ■ J ^Jp 2 + m 2 N f(p) d 3 p > M A , (10) 

i.e. the gas of nucleons has too large energy. This effect can be cured intro- 
ducing effective nucleon mass . 

A ■ J \jp 2 + (m^) 2 f(p) d 3 p = M A . (11) 

This corresponds to the reduction of the nucleon mass by 10-20%. 

Can the binding effect change our conclusions? The binding effect in 
nuclei is in general a complicated multiparticle problem and was solved in 
the exact way essentially only for three- and four-nucleon systems. The 
binding effect in heavy nuclei can be approximately included by introducing 

a notion of nucleon effective mass in the relation E N / A = \Jp 2 m/a + ( m N*) 2 - 
Since the binding effect is only a correction to our main effect of Fermi motion 
this approach seems justified. For illustration in FigJS]we show the result for 
nominal energy W = 20 GeV with nucleon mass modified by ± 150 MeV. 
This figure shows that the binding effect is not only qualitatively but also 
quantitatively a correction to our effect. This correction does not change the 
main trend of Raa{Pi)- A more refined calculation is not our intention and 
clearly goes beyond the scope of the present paper. 

How does the effect of Fermi motion depend on the type of particle mea- 
sured? In FigElwe present Raa for ix + + n~ (solid), K + + K~ (dashed) and 
p + p (dotted) at the center-of-mass energy W = 20 GeV and y — 0. In 
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this calculation Fermi gas model was used. The results for pions and kaons 
are similar, while the result for protons+antiprotons are somewhat different. 
The latter does not show the enhancement at large transverse momenta ob- 
served for pions and kaons. This is due to the parametrization for protons 
and antiprotons j2E] which is practically energy independent. Perhaps the 
difference for protons+antiprotons and pions/kaons can be used to identify 
experimentally the Fermi motion effect. 

In the present letter we have concentrated only on the effect of Fermi 
motion on nuclear modification factor. The Fermi motion in colliding nuclei 
is probably important in many other fixed-target (low energy) processes. For 
example global partonic QCD analyses use Drell-Yan data based on rather 
low-energy proton-nucleus collisions than elementary ones. The standard 
QCD approach is known to have problems in describing the transverse mo- 
mentum distribution of the lepton pair. In such analyses the Fermi motion 
of nucleons in nuclei is completely ignored. We expect that the effect of the 
Fermi motion should be dependent on the invariant mass of the two leptons. 
In general, the larger invariant mass, the larger Fermi motion effect may be 
expected. From the experience of the present paper one may expect large 
effects close to kinematical limits. Therefore even for the production of heavy 
gauge bosons (W,Z) in nucleus-nucleus collisions at RHIC one may expect 
sizeable effect due to the Fermi motion. These effects were never discussed 
in the literature and await a quantitative analysis. 

Finally we wish to mention another interesting, not fully explained, effect 
related to the CERES collaboration measurement of the elliptic flow param- 
eter V2 |2Z|- The existing hydrodynamical or thermodynamical models have 
serious problems in explaining a saturation of v% observed at larger transverse 
momenta (see e.g.|28j). It is quite probably that the Fermi motion effect may 
play important role in the full understanding of this effect. However, this ef- 
fect is very subtle and requires the inclusion of many effects in addition to 
the nucleon Fermi motion. This clearly goes beyond the scope of the present, 
intentionally simplified, analysis. 

4 Conclusions 

For the first time in the literature, in the present Letter we have analysed 
the role of the Fermi motion in nuclei in modyfing the so-called nuclear 
modification factor Raa- A large effect is found even at SPS energies. The 
details depend on the model of the Fermi motion as well as on the cross 
sections for elementary collisions. The latter are very difficult to calculate 
reliably, especially at low transverse momenta. Therefore we have made use 
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of a parametrization of existing experimental data. 

At SPS energies the Fermi motion causes a sizeable enhancement at p t > 
4 GeV. At RHIC energies the enhancement almost dissapears. 

This quantitative difference of the role of the Fermi motion at low and 
high energies is very important in the context of a quantitative understanding 
of jet quenching in dense nuclear media like quark-gluon plasma. At present 
the latter cannot be calculated from first principles and the model parameters 
are extracted by comparison to the experimental data. 

Even when alone, the Fermi motion leads to the Cronin effect. Different 
mechanisms like parton multiple scattering [Jj, initial transverse momentum 
broadening in nuclei [TJ |3] or perturbative saturation of gluon distributions 
jH] can contribute to the Cronin effect too. Our analysis shows that the Fermi 
motion effect is very important and must be included into the analysis before 
reliable conclusions about other nuclear effects may be drawn from such 
phenomenological studies. We find R AA (RHIC) < R AA (SPS) which makes 
the interpretation of the experimental data even more complicated. For 
example at the SPS energy W=17.4 GeV the inclusion of the Fermi motion 
would improve the agreement of the calculation in Ref. jl] with the WA98 
collaboration experimental data (see their Fig. 3). However, the inclusion 
of all effects like shadowing, quenching and Fermi motion in one consistent 
framework seems rather difficult at present. 

The Fermi motion effect seems to be one of a few different competing 
mechanisms modyfing Raa- Can one isolate experimentally the Fermi motion 
effect better? We think that the role of Fermi motion can be pinned down 
better in the collisions of light nuclei, where other effects are suppressed and 
the nucleon momentum distribution can be calculated from realistic models 
of nucleon-nucleon interactions. The examples are d — d or a — a scattering 
measured already in the past (see for example Better precision data 

at SPS energies would help to identify the effect of the Fermi motion via a 
comparison with realistic theoretical calculations 24 . Also very precise data 
on elementary collisions would be helpful when calculating the modification 
of Raa precisely. 
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Figure 1: Distribution of the true nucleon-nucleon cm. energies for Fermi 
gas model (dashed) and the Gaussian parametrization (solid) for nominal 
energy W = 20 GeV (left panel) and W = 50 GeV (right panel). 
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Figure 2: The energy dependence of the invariant elementary cross section 
for the reaction p + p — > ir° + X for different p t (=1,2,3,4,5,6,7,8 GeV) and 

y = o. 
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Figure 3: The nuclear modification factor R PhPh for W = 17.4 GeV and 
y = 0. The dashed line represent the result for the Fermi gas momentum 
distribution, while the solid line for the Gaussian parametrization of the data 
for the 12 C nucleus. The experimental data are from Ref.' L 8] for the 12.7 % 
most central collisions (N^n — 651). 
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Figure 4: Energy dependence of R\a at y = 0. The lines correspond to the 
nominal energies Wnn = 20, 50, 100, 200 GeV in the Fermi gas model (left 
panel) and using a simple Gaussian parametrization (right panel). The small 
fluctuations are due to poor numerics. 
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Figure 5: The effect of the effective nucleon mass variation for W=20 GeV. 
The solid mass is for standard free nucleon mass while the dashed lines are 
for nucleon mass modified by ± 150 MeV. In this calculation y = 0. 
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Figure 6: Fermi motion effect on nuclear modification factor for n + + 7r~ 
(solid), K + + K~ (dashed) and p + p (dotted) at W = 20 GeV and y = 0. 
In this calculation Fermi gas model was used. 
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